Templeting of Thin Films Induced by Dewetting on Patterned Surfaces 
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The instability, dynamics and morphological transitions of patterns in thin liquid films on periodic 
striped surfaces (consisting of alternating less and more wettable stripes) are investigated based 
on 3-D nonlinear simulations that account for the inter-site hydrodynamic and surface-energetic 
interactions. The film breakup is suppressed on some potentially destabilizing nonwettable sites 
when their spacing is below a characteristic lengthscale of the instability (A^), the upper bound for 
which is close to the spinodal lengthscale. The thin film pattern replicates the substrate surface 
energy pattern closely only when, (a) the periodicity of substrate pattern matches closely with the 
Xh, and (b) the stripe- width is within a range bounded by a lower critical length, below which no 
heterogeneous rupture occurs, and an upper transition length above which complex morphological 
features bearing little resemblance to the substrate pattern are formed. 

PACS numbers: 68.15.+e, 47.20.Ma, 47.54.+r, 68.08.De, 68.08.Bc 



Self-organization during dewetting of thin films on de- 
liberately tailored chemically heterogeneous substrates is 
of increasing promise for engineering of desired nano- and 
micro-patterns in thin films by templeting Jl], [| ||, ||, |[ ||, 
0, |[ ||, |l0[ [ll]]. On a chemically heterogeneous substrate, 
dewetting is driven by the spatial gradient of micro-scale 
wettability |Q , rather than by the non- wettability of the 
substrate itself. The latter occurs in the so called spin- 
odal dewetting on homogeneous sur faces H Hi- whilc 
the rupture of a thin film on a single heterogeneous patch 
is now well understood, patterned substrates pack a large 
density of surface features that are closely spaced. How 
does hydrodynamic interactions between the neighbor- 
ing heterogeneities affect the pattern evolution dynamics 
and morphology in thin films? This question, which is 
addressed here, is central to our understanding of how 
faithfully the substrate patterns are reproduced in a thin 
film spontaneously, i.e., how effective is the templeting 
of soft materials by dewetting route and what are the 
conditions for ideal templeting? An associated question 
for both the patterned and naturally occurring hetero- 
geneous surfaces is whether all the potentially dewetting 
sites remain active or "live" in producing rupture when 
they are in close proximity. These questions are resolved 
based on 3-D nonlinear simulations of the stability, dy- 
namics and morphology of thin films on periodic chemi- 
cally heterogeneous surfaces. 

The substrate considered consists of alternating less 
wettable and more wettable (or completely wettable) 
stripes that differ in their interactions with the overly- 
ing film. The key parameters of the substrate pattern 
are its periodicity interval (center-to-center distance be- 
tween two consecutive stripes, L p ) and the length-scale of 
the less wettable stripe (stripe- width, W) . The following 
nondimensional thin film equation governs the stability 
and spatio- temporal evolution of a thin film system sub- 
jected to the excess intermolecular interactions [ jl2| . 

dH/dT + V • [ff 3 V(V 2 ff )] - V • [# 3 V$] = (1) 



H(X,Y,T) is non-dimensional local film thickness scaled 
with the mean thickness h a ; $ = [2irh 2 /\A s \}[dAG/dH})-, 
AG is the excess intermolecular interaction energy per 
unit area and A s is the effective Hamaker constant for 
van der Waals interaction; X, Y are the non-dimensional 
coordinates in the plane of the substrate, scaled with 
a lengthscale (2TT^/\A s \) 1 / 2 h 2 ; and the non-dimensional 
time T is scaled with 127r 2 fJ.'jhg/ 'A 2 ; 7 and fi refer 
to the film surface tension and viscosity, respectively. 
The lengthscale of the spinodal instability on a uni- 
form surface is given by, A = (— iir 2 ^ / (d 2 AG / dh 2 )) 1 ^ 2 . 
On a chemically heterogeneous striped surface, $ = 
§(H, X, Y). At a constant film thickness, we model the 
variation of $ in X direction by a periodic step func- 
tion of periodicity, L p . Gradient of force V$, at the 
boundary of the stripes causes flow from the less wettable 
(higher pressure) stripes to the more wettable (lower 
pressure) stripes, even when the spinodal stability condi- 
tion 9$ / dH > is satisfied everywhere . It is known 
that a single stripe in the absence of its neighbors can 
cause rupture only if its width exceeds a critical length- 
scale, W c « A |J. 

We consider general representation of antagonistic (at- 
tractive/repulsive) long range and short range inter- 
molecular interactions applicable to the aqueous films 
(eq. 2) ^] and polymer films on high energy sur- 
faces like sillicon (eq. 3) j|, |l2], |l6). In both the cases 
the long range van der Waals force is stabilizing and thus 
films thinner (thicker) than a critical thickness are unsta- 
ble (metastable) . 

AG = -(AslYl-nh 2 ) + S p exp(-h/l p ) (2) 

This potential represents aqueous films when the effective 
Hamaker constant, A s and S p are both negative (long- 
range van der Waals repulsion combined with shorter- 
range attraction, e.g. hydrophobic attraction) JlTj . On 
the hydrophobic stripes (A* = -1.41 x 1CT 20 J, S* = 
—65 mJ/771 2 and l p = 0.6 ran), the film is spinodally 
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FIG. 1: Equilibrium thin film profile on a striped surface 
(W = 0.16A). The less-wettable stripes are denoted by black 
rectangles. Figs la, lb and lc is for a 6nm film for increasing 
periodicity. Fig Id shows the profile for a 5.5 nm thin film on 
a same substrate shown in Fig. lb. 



unstable (d$/dH < 0) below a critical thickness, h c of 
7.4 nm. The equilibrium contact angle obtained from 
extended Young-Dupre equation, cosO = 1 + AG(/i e )/7 
is 61°; where h e is the equilibrium thickness. On the 
hydrophilic stripes {A s = -1.41 x 10~ 20 J, S p = -0.61 
mJ/ra 2 and l p = 0.6 nm), the film of any thickness is 
spinodally stable (d$/dH > 0) and completely wets the 
surface. 

Equation (1) was numerically solved using a central dif- 
ference scheme in space combined with Gears algorithm 
for stiff equations for time marching and periodic bound- 
ary condition. The critical length, W c of a single hetero- 
geneity (hydrophobic patch) that engenders rupture on 
a large hydrophilic substrate in this case is 0.08A. 

Ideal templating requires that dewetting should occur 
on every less wettable site on the patterned substrate. 
Fig. la shows a single rupture on a substrate contain- 
ing four potentially dewetting sites (ns — 4) when the 
periodicity interval is small, L p = 0.45A. The result- 
ing lone drop spans across the remaining "inactive" hy- 
drophobic sites. Fig la to lc show increased number 
of dewetted stripes {no — 1, 2, 4) and resulting liq- 
uid domains cis Lp is increased from 0.45A — > 0.68A — > 
0.9A. This, as well as a large number of simulations 
not displayed here showed that on a substrate contain- 
ing many potentially destabilizing sites, only the sites 
(randomly picked) separated by a characteristic length- 
scale of the instability, Xh, of the order of the spinodal 
scale, A remain "live" or effective in causing the film 
rupture. Dewetting on the remaining intervening sites 
is suppressed since rupture on each site would require 
surface deformations on smaller scales resulting in high 
surface energy penalty. The ratio of this characteristic 
lcngthscale, Xh, to the spinodal length scale, A, obtained 
from simulations, decreases as the potential difference 
across the stripe bounadry, (A$) increases (Fig. 2), 
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FIG. 2: Variation of the ratio of characteristic to spin- 
odal length scale with potential difference A<6, introduced by 
changing the S p of the more wettable stripes. —d<f>/dh, for 



curves 1 (h = 6nm, S p 



-65mJ/m ), 2 (h = 4.5nm, 



Sp = -19.5mJ/m 2 ) and 3 (h a = 4.5nm, Sp = -39mJ/m 2 ) 
are 6.47xl0 12 , 2.4xl0 13 and 5.4xl0 13 J/ra 4 , respectively. 

where A$ = (& h — <f>), evaluated at the initial thick- 
enss. For a large value of spinodal parameter (e.g. for 
thinner films), the ratio approaches very close to 1 unless 
A$ is very large (curves 2 and 3 of Fig. 2). A simple 
scale analysis (H = 1 + e, V 2 <I> ~ A$/cL 2 , c < 1), of Eq. 
(1) 1 12 leads to the following characteristic lengthscale 
for growth of instability: X h cx [A$/ec — d$/dH]~ 5 . 
Thus the ratio of heterogeneous to spinodal lcngthscale, 
X h /X cx (1 + A$/(-ec<9$/<9ff))- a5 ; and X h ~ A when 
spinodal term —dQ/dH is very strong compared to the 
applied potential difference. For very large potential 
difference, Xh/X cx (AcA) -0 ' 5 . Thus the upper limit of 
the characteristic lengthscale of the instability is always 
bound by the spinodal length, A. An interesting impli- 
cation is that even when the rupture occurs by the het- 
erogeneous mechanism on a substrate containing a large 
density of heterogeneities, the lengthscale of the result- 
ing pattern can give the illusion of spinodal dewetting 
by mimicking the characteristic lengthscale of the lat- 
ter. Differentiating true spinodal dewetting from het- 
erogeneous dewetting therefore requires a careful consid- 
eration of their distinct time scales and morphological 
features |12||. Further, the contact line always remains 
pinned at the boundary of the stripes whenever alternate 
stripes are completely wettable (Figures la — Id). 

Decreased mean film thickness increases both the 
terms, A$ and —d&/dh and thus decreases Xh nonlin- 
eraly. Thus, a greater number of nonwettable hetero- 
geneities become active in causing the film breakup for 
thinner films. Comparison of figs, lb (thickness = 6 nm, 
A = 0.67 /zm) and Id (thickness = 5.5 nm, A = 0.42 /jm) 
clearly shows this for an identical substrate pattern. This 
transition uncovered based on dynamical simulations is 
also in conformity with the earlier equilibrium energy 
considerations |] that are related however only to the 
equilibrium structures. 

From the above discussion, a necessary condition for 
good templating is that L p ~ Xh(< A) which ensures 
dewetting on every less wettable site (ud — n>s)- The 
condition L p ~ A^, for njj = ns, also remains valid 
when both the stripes are nonwettable. However, in this 
case the contact line can move across the boundary of 
the stripes to the more wettable part. We simulated 
the 3D morphologies during the evolution of a polymer 
like film on an oxide (low-energy) covered silicon (high- 
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FIG. 3: Morphological evolution in a 5 nm thick film on a 
striped surface (W = 0.8 /im = 0.53 A) In Figures 3A to 3D 
periodicity of surface pattern, L p are 1 fim, 1.2 /im, 1.45 /im 
and 3 /im respectively. The first image in this figure as well 
as in the subsequent figures represents the substrate surface 
energy pattern; black and white represent the more wettable 
part and less wettable part respectively. For other images 
describing the film morphology, a continuous linear gray scale 
between the minimum and the maximum thickness in each 
picture has been used. 



energy) substrate. An analytical representation of com- 
bined antagonistic (attractive/repulsive) long and short- 
range intcrmolecular interactions for a polymer-like film 
on a coated (e.g. oxide covered) substrate isjx3] : 



12ttAG = [(A a - Acx)/(h + d cl + d c2 f 
+ (A cl - A c2 )/(h + d c2 f + A c2 /h 2 ] 



(3) 



Negative value of the effective Hamaker constant on 
the substrate [A s = A% = -1.88 x 1CT 20 J] signifies a 
long-range repulsion, whereas a positive value on coating 
[A c i = A ] l x = 1.13 x 10~ 20 J] represents an intermediate- 
range attraction |T3|. The non- wettable coating (e.g. ox- 
ide) thickness [d c \ = 2.5 nm, = 4 nm] is increased 
on alternating stripes which causes the macroscopic con- 
tact angle to increase from 0.58° on the more wettable 
stripes to 1.7° on the less wettable stripes. A still shorter- 
range repulsion may arise due to chemically adsorbed 
or grafted layer of the polymer (A C 2 = —0.188 x 10~ 20 
J, d C 2 = lnm are the film Hamaker constant and the 
thickness of the adsorbed layer, respectively)) 12 . The 
representation, eq.[3], is chosen for illustration without 
affecting the underlying physics, which we have verified 
for many different potentials. 

Figures 3A — 3D depict the effect of periodicity in- 
terval, L p on the self- organization of a 5 nm thin film 




FIG. 4: 3D morphologies at late stages of evolution for films 
of thickness (A) 2.5 nm, (B) 6.9 nm and (C) 8 nm on a striped 
surface with W = 1.2 /im and L p = 3 /im. 



that is spinodally unstable on both type of stripes. In 
each case, the evolution starts with local depressions on 
the less- wettable stripes (Figs. 3A — 3D). For L p suffi- 
ciently smaller than X(L p = 1/im, A = 1.51/im based on 
7 = 38 mJ/m 2 and /j, = lkg/m.s), the isolated holes or 
depressions that form along the less wettable stripes grow 
onto the more wettable regions and coalesce with each 
other rapidly (Fig. 3A) leading to a disordered struc- 
ture and very poor templeting. Increasing periodicity 
(L p = 1.2/xm; Fig. 3B) leads to more ordered dewetting, 
but the number of dewetted regions remain less than 
the number of less wettable stripes, and defects evolve 
at late times (e.g., holes in the liquid ridge; image 4 of 
fig. 3B). For L p > A (Figures 3C and 3D), the num- 
ber of dewetted stripes equal the number of less wettable 
stripes. However, templeting in the form of liquid ridges 
with straight edges is best at an intermediate periodicity, 
L p ~ Afr(< A) (Fig. 3C). A further increase in L p makes 
the width of dewetted region bigger than the width of 
the less wettable stripe (image 4 of Fig. 3D), i.e., the 
contact line resides in the interior of the more wettable 
stripe, rather than close to the boundary. 

The transition of surface morphology can also be 
clearly understood from Fig. 4, which shows that ideal 
templeting occurs for an intermediate thickness (3rd im- 
age of Fig. 4) for which the spinodal wavelength, A(= 
3/im) is equal to L p . Thinner films evolve by the forma- 
tion of droplets on the more wettable stripes due to the 
spinodal mechanism assisted by the Rayleigh instability. 
For thicker films (A > L p ), dewetting occurs on fewer 
stripes with the formation of holes in broad liquid ridges 
(e.g., last image of Fig. 4). 

We found that besides L p , the other parameter that 
governs the thin film pattern is the stripe- width (W). 
It is already known that for a single heterogeneity on a 
large substrate, increase in the width of the heterogene- 
ity shifts the onset of dewetting from the center of the 
heterogeneity to the boundary of the heterogeneity |L2| . 
Figs. 5A — 5D shows the transition of patterning with 
the increase of the stripe- width (W), keeping the peri- 
odicity, L p fixed (> A). For small stripe-width (~ 0.4A; 
Fig. 5A), rupture is initiated at the center of less wettable 
stripes by the formation of depressions that coalesce to 
form rectangular dewetted regions of widths greater than 
the stripe- width. For a larger stripe- width (W ~ 0.8A), 
dewetting is initiated by a layer of holes at each of the 
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FIG. 5: Morphological evolution in a 5 nm thick film on a 
striped surface (L p = 3 fim ~ 2A). W = 0.6 /Ltm, 1.2 fi m, 2.1 
fim and 2.7 /im, respectively for figures 5A — 5D. 



two boundaries of the stripe (image 2 of Fig. 5B). Coa- 
lescence of these two layers of holes leads to the dewet- 
ted regions containing some residual droplets (image 4 of 
Fig. 5B). Further increase in the stripe width leads to 
the formation of two layers of holes on each stripe sepa- 
rated by an elevated liquid cylinder (image 2 of Fig.5C). 
Thus, at an intermediate stage of evolution, the number 
of cylindrical liquid ridges becomes twice of the number 
of more wettable stripes on the substrate (image 3 of 
Fig. 5C). Eventually, liquid ridges on the less wettable 
region disintegrate into droplets (not shown). Laplace 
pressure gradients cause ripening of the structure lead- 
ing to the merging of droplets with the liquid ridges (im- 
age 4 of Fig. 5C). Further increase in the stripe- width 
breaks the whole order of the substrate pattern due to 
the formation and repeated coalescence of several layers 
of holes on each stripe that eventually evolve into arrays 
of irregular droplets. Thus, a good synchronization of 
the thin film morphology with the substrate pattern re- 
quires an upper limit on the stripe-width (W* ~< 0.8A, 
for this case). Based on a large number of simulations 
for different film thickness and for different systems, we 
have verified that a breakdown of templeting due to the 
off-center dewetting occurs at stripe- widths in excess of 
about 0.7A - 0.8A. 

In conclusion, the 3-D thin film morphology on a pat- 
terned substrate during dewetting can be modulated pro- 
foundly by a competition among the time scales (spinodal 
and heterogeneous) and length scales (spinodal, stripe- 
width, periodicity, thickness) of the problem. On a 
striped substrate, ideal templeting in the form of cylin- 



drical liquid ridges covering the more wettable stripes 
occurs when (e.g., image 4 in fig. 3C): (a) periodicity 
of the substrate pattern is very close to the character- 
istic lengthscale, Xh , and (b) the stripe width is larger 
than a critical width, W c , which is effective in causing 
rupture by the heterogeneous mechanism, but smaller 
than a transition width (W t ~ 0.7A — 0.8A), that en- 
sures initiation of dewetting at the stripe- center. Even a 
very weak wettability contrast (e.g. Fig. 3) is sufficient to 
align the thin film pattern with the templete under the 
above conditions. Predictions of our simulations show 
close resemblance to the recently reported experimental 
observations on dewetting of polymers on patterned sur- 
faces that reveal: (a) best templeting for an intermediate 
thickness film[|| and (b) correspondence between the 
natural length scale and periodicity of the substrate pat- 
tern for good templeting H, |], . Finally, on a completely 
wettable substrate containing many potentially destabi- 
lizing sites, only the sites (randomly picked) separated 
by a characteristic lengthscale (of the order of the spin- 
odal scale), Ah remain "live" or effective in causing the 
film rupture. It is hoped that this study will help in the 
design and interpretation, creation and rational manip- 
ulation of self organized microstructures in thin films by 
templeting. 

Discussions with G. Reiter are gratefully acknowl- 
edged. 
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